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Abstract 

We calculate the perturbative correction to every cluster coefficient of a gas of anyons through 
second order in the anyon coupling constant, as described by Chern-Simons field theory. 
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It has been known for some years that two dimensional physics offers the possibihty of 
fractional statistics (anyons) interpolating between bosons and fermions [1]. Either for its 
possible relevance in some areas of condensed matter physics, or for the deep implications 
it has on the concept of identical particles, the study of the thermodynamical properties of 
the anyonic system has deserved great effort since the past decade. However, until recently, 
the low density, high temperature expansion could only provide the second virial coefficient 
[2]. Some restricted information about the third virial coefficient was obtained in refs. [3,4]. 
Numerical computations have also been performed [5-7]. The root of the problem is the 
unsolvability of the A'"-anyon quantum mechanical system beyond N = 2. 

This obstacle can be overcome making use of second quantized theory. The price one 
has to pay is that perturbative field theory only provides an expansion in powers of 9 (the 
statistical parameter) near bosonic or fermionic statistics. 

In this approach, a problem is posed by the singular nature of the short range statistical 
interaction. In refs. [8-11] a solution is proposed introducing a two-body non-hermitian 
interaction. Divergences must be carefully controlled with harmonic regulators or boundary 
conditions. In this way, Dasnieres de Veigy and Ouvry [11] have been able to compute the 
partition function to order 9^, expanded to all orders in the fugacity z — exp{fj,/T). 

Recently another model has been considered to solve the problem [12]. There, fractional 
statistics is implemented in standard way through non-relativistic matter interacting with a 
Chern-Simons gauge field. UV divergences are cancelled with a two-body 5-function potential 
when the starting statistics is bosonic. Starting from fermionic statistics does not require 
such interaction. This model is used for to the computation of the second virial coefficient, 
to three-loop order. In this letter we develop further the model and obtain the solution to 
all orders in z, to the same loop order. Remarkably, our results are in full agreement with 
those of refs. [10,11] t. 

^ The relation between both approaches is discussed in ref. [12]. 
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The setting of our model, as proposed in ref. [12] consists of a spinless non-relativistic 
self-coupled matter field interacting with Chern-Simons gauge bosons. The Lagrangian 
density will be given by 

111 1 

£(t,x) = -—dtaixa + -aoB - — (V • a)^ +^UDoij-— |D^|2 +;,^t^ _ ^ (^t^)2^ 

Do^dt + iao, D = V-ia. (1) 

Here 5 = V x a is the magnetic field, n is the chemical potential and p is a gauge fixing 
parameter. The Coulomb gauge to be used refers to the choice p — 0. The anyon coupling 
K, can be related to the statistical angle parameter 9 as k, — 29. Also, the strength of the 
contact interaction a must be taken [12] as 

a = {l + C)-, (2) 
m 

{( = 1 for bosons; ( = —1 for fermions). With this choice, the divergences appearing 
for fiducial bosonic statistics are exactly cancelled; for fermions, a vanishes, not needing 
additional interactions. 

In the imaginary-time formalism of thermal field theory [13], the functional integral 
expression for the grand partition function involves an integration over imaginary time from 
to /3 = , 

Z(ii,T,V)^ f VaoVaV^^^V^l; ex.p( f"^ dr f d'^xC(t^-iT,:x.)\ . (3) 

J (anti)periodic ^Jo J ' 

Here (anti)periodic means that the integration over fields is constrained so that ■?/'(x, /?) = 
±'0(x, 0) where the (lower) upper sign refers to (fermions) bosons. 

Now we can proceed to expand in a power series in k, by using a set of diagrammatic 
rules which follows from (1). These are listed in Table 1. The diagrams describing the 
perturbativc scries for InZ have the form of connected closed loops. It should be noted that 
there will be a factor of ^iA left over for each graph, corresponding to the extensivity of In Z. 
This factor cancels out in expressions for the pressure. 
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The contribution to In Z up to order involves a set of three-loop diagrams. The non- 
zero graphs are shown in fig.l ^ Note that all graphs with only one gauge boson line vanish 
due to the index summation associated with vertices. In order to keep the correct order of 
the operators according to their r-values, one must insert a factor e*'^"^ whenever a particle 
line either closes on itself or is joined by the same instantaneous interaction line. We take 
?7 — > 0+ at the end of the calculation. This implies the following basic frequency sums 



-C ,„„ > ^-: = -Cnq, (4) 



/J^iu-„ + Ai-qV(2m) %xp [,3 (qV(2m) - /I.) ] -C 

-y 



/3^za;, + /x-q2/(2m) 2 
As a consequence, graphs (c), (g) and (h) are of order . 

The above explained formalism allows us to derive the perturbative corrections to the 
pressure 

P = ^ In Z. (6) 

To each order in k our results will give the pressure as a function of the fugacity z — ex^{^/T). 
This corresponds to the cluster expansion 



^ oo 



(At = (27r/mT)^/2 is the thermal wavelength). However, the statistical dependence is more 
clearly manifest in the virial expansion 

^^-^f^a,{nXl)\ (8) 

with 

„.iVM = 4@. (9) 

Both expansions can be related easily order by order: 

ai = 1 , 

2 ^ ^3 



a3 = 4a^-2^, (10) 



^ Graph (h) was missing in ref.[12]; since its contribution is of order the results of ref.[12] remain 

unaffected. 
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We proceed now to calculate the first two orders of the perturbative expansion in n. 



First order 

As argued before, to this order, diagrams with one gauge boson line vanish. This leaves 
only diagram (a). As a result, fermionic statistics receives no corrections of order k. The 
contribution to the pressure can be written 

'2.^ . n 2 



-<"-f[ii:/^av,p)e 



-~ 111^(1 -^). (11) 



This result agrees with that of refs. [9,11]. It is easy to show from eqs. (10) and (11) that 
only the second virial coefficient is corrected to this order. 

Second order 

Diagrams (b), (e) and (f) are hardest to compute. They are given by 

.2 r ^2„ j2i, j2 



L a f d p d k d a ^r., . ^n, 

^' = 8^ Y. J p) «°(-' + p + 

X g^{uj2,k) ^°(a;2 + i^,k + q) 



/t^ f d^p d^k d^q ripTik+q + ^Cn^ripni^+q 
' ^'2mj (27r)2 (27r)2 (27r)2 q • (p - k) ' ^ ' 



T^p > 1^ f d p d'^k rf^q , ^n, , 

^ (2JF ^ ^ + p + 1' 

2 

xg\u;2,k)g\u;2 + u,k + ci) " ^ ^] 



(P - k)2 q2 

/• d^p d^k d^q ripTik+q + 2Cnknpnk+q [(p - k) x q] 
^ 2m i (27r)2 (27r)2 (27r)2 q • (p - k) (p - k)2 q2 ' 



Or. , -.1. , „^ iP-k) -q 



(P - k)2 q2 

Prom these equations it is straightforward to notice the following relation 

pb^-{l + C){P' + pf), (15) 

reflecting the cancellation of exchange diagrams in the bosonic case. Identity (15) implies a 
substantial reduction of calculations; to find + P'^ + , the only integral to compute is 
(12). 

There are two different terms in (12): one involving the product of two n's and another 
with three n's. They give quite distinct contributions: 

2 rx* 
pb^pe^ pf = ^ ln2(l _ ^z) / dx $(x) 

47rmAy Jq 

Here x* is a UV cutoff and $(a;) is the plasma dispersion function [14], 

^x) = 2 dt^=== = 2e-- / dte' , (17) 

Jo y X V J 

its derivative being given by 

^'{x) = 2- 2x^x) . (18) 

The UV divergence comes from the loop with the seagull vertex and from the divergence 
in the second Born approximation for scattering with a (5-function interaction [15]. These 
divergences will cancel when we consider the whole set of diagrams. 

We evaluate now the ring contribution from graph (d) of fig. 1, 

ud_^^^ f d\ n%o(;/,q)S°(i/,q) 
^-2^^7(2^ ? ■ 
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where Il^oo{u,q) and S°(i/, q) represent the density-density correlator and the transverse 
component of the current-current correlator respectively at the lowest order, 

n%-(^,q) = -;^^E /(^[(P + i)^(P+f).-^Vi,P)^Vi + ^,q + p) 

= -n%o (-n, q) I ^ + S° (.n, q) {SiJ - ^) ■ (20) 

Performing the sums and integrals, it yields (the first part comes from two-n terms; the 
second, from three-n terms) 

.2 rx* 



t tj Th — -L S — J- 

The ring loop with a seagull vertex has yielded a UV logarithmic divergence. This exactly 
cancels the infinities we have found before. To see how this happens, we quote some results 
about the plasma dispersion function, namely 

r dx e-'*^' - - + 2(1 + a)$(a;)) = 2 , (22) 

Jo \ X X / 



Hdx e-^ ^x) = -^ln(^^±l) . (23) 

Both of them follow from eqs. (17) and (18). Now, first use eq. (22) (with a = 0) to compute 
the two-n contribution from eqs. (16) and (19): 



K 



2 



ln'(l - C^) ■ (24) 



Split the three-n term in eq. (16) in two parts: one is what must be added to the three-n 
term in eq. (21) so that eq. (22) can be used. The sum in z can be performed: 



ln(l - C^) (^^T; + 1^(1 - C^)) • (25) 



The remaining part of eq. (16) is integrated using eq. (23): 



2 1—2 l~n—l „ , 

r^E(c-)T E f^in(§^), (26) 



where 



_ n(2g-/ + n) ^ / (/-8)(n + 5) 

s2(/-s)(i-n-s) ' ^^^^-y s(i-n-s) ' ^ ' 

Due to the logarithms, it does not seem possible to perform the sums. Notice also that these 

terms give no correction to the first three cluster (and virial) coefficients. 

It only remains to evaluate diagrams (c), (g), (h). These are the easiest to compute. 
They yield finite results: 

2 

^TT^T-^ 1^2(1 - C^) , (29) 



P 



h 



The solution we have found consists of a summable part (eqs. (24), (25), (28), (29) and 
(30)) and a non-summable logarithmic part, eq. (26). Collecting the summable terms we 
get 

We find complete agreement with the results of ref. [11]. The third virial coefficient obtained 
from these results does not depend on C; hence it satisfies the 'mirror symmetry' discovered 
by Sen [3]. For higher virial coefficients, this symmetry docs not hold. 

To conclude, we have shown that the model based on standard perturbative Chern- 
Simons field theory coupled to non-relativistic matter accounts for the thermodynamic prop- 
erties of the anyon gas. As an illustration, we have derived a closed formula for the n-cluster 
coefficient through second order in the anyon coupling constant. 
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Figure Captions 

Fig. 1. Non zero diagrams contributing to the pressure to the second order in n. The sign ± 
refers to Bose or Fermi propagators. Graphs (c), (g) and (h) do not contribute to the second 
virial coefficient. Combinatoric factors are shown in the diagram. 
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